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Abstract 

We explore the low-temperature behavior of the Abelian Higgs model in AdS 4 , away 
from the probe limit in which back-reaction of matter fields on the metric can be neglected. 
Over a significant range of charges for the complex scalar, we observe a second order phase 
transition at finite temperature. The symmetry-breaking states are superconducting black 
holes. At least when the charge of the scalar is not too small, we observe at low temperatures 
the emergence of a domain wall structure characterized by a definite index of refraction. We 
also compute the conductivity as a function of frequency. 
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1 Introduction and summary 

The Abelian Higgs model in four-dimensional anti-de Sitter space (AdS^) is specified by the 
action 

R-\Fl u -\{d»-iqAM 2 -VM) , (1) 

where V depends only on the magnitude of the complex scalar field ip, not its phase. Among 
the simplest solutions to the classical equations following from the action (jl]) is the AdS$- 
Reissner-Nordstrom black hole solution, hereafter RNAdS: 

ds 2 = e 2A {-hdt 2 + (dx 1 ) 2 + (dx 2 ) 2 ) + ^ 

A = — h = l- eLn 2 e-^ L + ^ e ~^ L ^ 
L 4 

$ = p K 2 (e- rH/L - e~ r/L ) ij) = . 

In pQ, following earlier work [2], it was suggested that there are black hole solutions 
to the classical equations of motion following from (jl]) which spontaneously break the U(l) 
gauge symmetry associated with phase rotations of ip. The symmetry-breaking solutions can 
be thought of as superconducting black holes. More precisely, there is a superconducting 
condensate of the scalar field ip which floats above the horizon. A calculation of the critical 
temperature T c below which black hole superconductivity occurs was outlined in [2] and 
carried through for a few choices of parameters such as q and the anti-de Sitter space radius 
L. The conformal field theory (CFT) dual to a superconducting black hole of this type 
is a symmetry-breaking phase whose order parameter is (O^), where is the operator 
dual to ip. In the simplest setup, the phase rotations of are a global symmetry of 
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the boundary theory; however, one can weakly gauge this U(l) symmetry, and then the 
symmetry-breaking in the boundary theory can also be described as superconductivity. This 
transition from breaking a global symmetry to breaking a gauge symmetry has a precise 
parallel in BCS theory: all the dynamics of electrons and phonons can be carried out in 
ignorance of the gauge interactions; but once the condensate Cooper pairs are formed, one 
can go back and note that it breaks the abelian gauge symmetry. As has been emphasized in 
[3] , many — though not all — of the macroscopic features of superconductors can be understood 
as a consequence of an effective field theory treatment of the breaking of the abelian gauge 
symmetry. 

In |3], the picture of superconducting black holes was fleshed out in the limit of large q, 
where the matter fields do not back-react on the geometry. It was further shown in [3] that 
for the choice 

V = -^ + m 2 \^\ 2 , (3) 

with m 2 L 2 = —2, the conductivity defined from a two-point function of the conserved current 
dual to exhibits gap-like behavior. This feature, remarkable for its qualitative (and 
perhaps semi-quantitative) similarity to real-world superconductors, cannot be anticipated 
on general grounds of symmetry-breaking and effective field theory. However, as observed in 
[3], the large-g limit doesn't commute with the small-temperature limit. So an outstanding 
question is what the ground state of the system is. 
In [5], a ground state was proposed for the choice 

V- = -^+m 2 h/f + ?^| 4 , (4) 

with m 2 < and u > 0. On the gravity side, this ground state has the structure of a 
domain wall. On the ultraviolet side is AdS± with ip = 0. On the infrared side is Ad S4 
with \ip\ = \J—m?/u. In between is a charged condensate of ip. Asymptotically far into the 
infrared, there is no electric field. In terms of the dual field theory, this means that all the 
charge is carried by the condensate: none remains in the normal state. This is a satisfying 
picture, but there are some strange features. The infrared AdS& signals the emergence of 
relativistic conformal symmetry, SO (3, 2), in the infrared. One consequence of this symmetry 
is that there is a maximum speed of propagation for signals at far infrared energies, and this 
speed is less than the speed of light. In other words, the geometry has a non-trivial index of 
refraction for infrared signals. A second consequence is that the real part of the conductivity 
exhibits power-law decay for small u: Re a oc u s with an exponent 5 that varies continuously 
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as a function of q, m 2 , u, and L. Power-law decay at small frequency means that there is no 
gap. It was suggested in [5] that for a more general class of potentials, a similar domain wall 
structure would arise, but the infrared side would not be anti-de Sitter. Instead it would 
be a warped geometry whose isometry group is the Lorentz group SO(2, 1) (together with 
translation invariance in space and time). This would again imply a non-trivial index of 
refraction, but the conductivity would not be forced to exhibit power-law scaling at small uj: 
it might instead exhibit gapped behavior. 

The aim of the current paper is to go beyond the linearized calculations of [1] and the 
probe approximation of [J] and track the thermodynamically preferred phase of the Abelian 
Higgs model, from T c down to low temperatures. This task proceeds largely by finding 
numerical solutions to the equations of motion, of the form 

fa* = e 2A(r) [ h{r)dt 2 + (dx l)2 + ,^2)21 + ^_ 

1 h(r) (5) 

$ = $(r) ip = tp(r) , 

where $ = A$ is the timelike component of the gauge potential. We always insist that 
h(r) should have a simple zero at a finite value of r: that is, there is always a regular, 
finite-temperature horizon in our solutions. In order for A^dx^ = $>dt to be well-defined at 
the horizon, we must have $ = at the horizon. We also insist that ip oc t~^ A near the 
boundary of AdS^, where A^ is the larger root of 

A^(A^ - 3) =m 2 L 2 . (6) 

This condition is equivalent to specifying that the lagrangian of the dual field theory isn't 
deformed by C^Q 

We have two main findings: 

1. Superconducting black holes of the form ^ are thermodynamically preferred over 
RNAdS black holes below a temperature T c at which a second order phase transition 
occurs. This was speculated to be the case in pQ and asserted to be true in the probe 
approximation in [I] ; also there was a claim in [7] of unpublished calculations verifying 
that superconducting black holes are favored. More recently, the phase transition was 



1 For a range of dimensions, 3/2 < A^, < 5/2, it is consistent to choose instead the smaller root of pj. 
This choice corresponds to a different CFT, whose correlators can be systematically related to the original 
CFT [6]. A discussion of superconducting black holes with these alternative boundary conditions has been 
given, in the probe approximation, in [4]. 
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studied in the probe approximation in the more general setting of non-zero superfluid 
velocity, and a tricritical point was found [HI E]. While this paper was nearing com- 
pletion, we received [10J, which has some results overlapping ours on thermodynamics 
away from the probe approximation. The results reported here together with the lit- 
erature just cited present quite a convincing case that superconducting black holes 
describe the preferred phase below the temperature T c obtained from analyzing linear 
perturbations around RNAdS. 

2. At least for large enough values of qL, an 5*0(2, 1) symmetry and a non-trivial index 
of refraction arise at sufficiently small temperatures. This is associated with charge 
being completely expelled from the black hole, or in field theory terms, to charge 
being completely carried by the condensate. As qL decreases, the index of refraction 
decreases. For sufficiently small qL, it is difficult to say whether S0(2, 1) symmetry 
arises in the limit of small temperatures. If it does, the relevant temperature scale is 
very small indeed. 

We have been unable to construct symmetry-breaking solutions at zero temperature. Such 
solutions would be dual to the ground state of the conformal field theory in the presence of 
a finite charge density. It is reasonable to expect that such solutions exist, and that they 
exhibit an emergent Lorentz symmetry in the infrared, at least when qL is not too small. 

The rest of this paper is organized as follows. In section [2] we collect analytical results 
on the equations that determine the background as well as equations that can be solved to 
find the conductivity. In section [3] we report on the result of our numerical studies of these 
equations, focusing on thermodynamics, a finite index of refraction, and the conductivity. 



2 Analytical results 



Plugging the ansatz (|5| into the equations of motion resulting from ([I]), one finds the following 
second order differential equations: 




ti' + 3A'ti = e~ 2A & 2 + V (8) 

he 2A 

q>» + A'& = ^-$^ 2 (9) 



(10) 
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where primes denote d/dr. There is also a first order constraint, which if satisfied at one 
value of r must hold everywhere, provided the equations of motion ([7 10) are also satisfied: 



h^'i + e - 2A q 2 ^ 2 - -he- 2A & 2 - 2hh'A' - Qh 2 A' 2 - hVty) = . (11) 

We assume that ip is everywhere real. This makes sense because it costs energy for its phase 
to vary, and with variation of fields only in the r directions, there is no non-trivial topology 
for its phase to wind around. 

Precisely the equations (|7|)-([TT| were derived in [3j. But unlike in that paper, here we 
require that h has a simple zero at a finite value of r. Using the freedom to shift r by an 
additive constant, we may require that the horizon occurs at r = 0. Using also the freedom 
to rescale t and x = (x 1 ,^ 2 ) by separate multiplicative factors, we may also require 

A{0) = and h'{0) = 1 . (12) 

In order to specify a solution to the equations of motion with a regular horizon, we need 
only two further conditions: 

ip(0) = and $'(0) = $i . (13) 

A standard parameter-counting argument suffices to show that there is now a well-defined 
Cauchy problem. Here is how the argument goes. The four second order equations have 



eight integration constants, but one is used up by the constraint (11), one more is used up 



by insisting that the horizon is at r = 0, and four more are used up by the explicit conditions 



( 12 )-( 13 ). This appears to leave two constants undetermined. But in fact, the existence of 
a regular horizon implicitly requires that $ = at the horizon and that ip is regular there, 
and these two conditions amount to fixing the final two integration constants. So, based on 
an analysis of the horizon boundary conditions, there is a two-parameter family of solutions, 
parameterized by an d $i- 

There is one additional boundary condition at infinity: as we have already mentioned, 
we require that ip cx e~ A ^ A near the boundary of AdS±. This requirement amounts to a non- 
linear constraint on ($ 1; if) ). Any pair of values that satisfies the constraint corresponds to a 
black hole that spontaneously breaks the abelian gauge symmetry: that is, a superconducting 
black hole. In practice, one must make some further restrictions on the allowed class of 
solutions: for large r, A must asymptote to a linear function of r with positive slope, while 
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h and $ must asymptote to constants. There are solutions which violate one or more of 
these restrictions: for example, A — ► — oo as h — > +00 is fairly common. We believe all such 
solutions are singular. Certainly they are not asymptotically AdS^, and as such they can 
naturally be excluded. By studying series expansions near the boundary of solutions that 
are asymptotically AdS^, one can show that 

A(r) = a x r + a + . . . 
h(r) =H + H 3 e- 3A + ... 

(14) 

$(r) =p +Pie A + ... 
ip(r) = e _A * A + . . • , 

where in each equation, . . . indicates terms which are subleading at large r to the ones shown. 
Because r = is by assumption the horizon, i.e. the largest value of r where h has a zero, 
we must have H > 0; also, as remarked previously, we require a\ > 0; but all the other 



constants in ( 14 ) might in principle have either sign. In practice, we generally find po > 
and pi < 0. 

We make one further restriction on the solutions we study below: ip(r) is not allowed to 
have any nodes. This is based on the idea that oscillations in ip can only add to the energy, 
and there is no topology to support them, so a solution with nodes is probably unstable 
toward decay to a solution with no nodes. Such expectations have been borne out in related 
studies of solutions of the Einstein- Yang-Mills equations: for a recent review see [TT] . 

Having found a solution to ([7|-( 11 ), one can determine its thermodynamics by extracting 



the coefficients shown in (14). The energy density, entropy density, temperature, chemical 



potential, and charge density can be read off as 



g 3 c = ^ e 2A(0) 



k 2 LHq k 2 



T=—e A ^^=L (15) 



Pi f rr, 

P = wr^ f = e-Ts. 



We adhere to conventions of [7], which discusses the thermodynamics of a model similar 
to the Abelian Higgs in AdS±. Energy density, temperature, and chemical potential are 
measured with respect to the Killing time y/H^t instead of t. We also follow [7] in defining 
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rescaled thermodynamic quantities with no k dependence, as follows: 



n 2 e A k 2 s 



(2tt) 3 L 2 (2tt) 3 L 2 

(16) 

» K P i « / 

P (2tt) 3 L 2 J (2tt) 3 L 2 ' 

The states in the dual conformal field theory that charged black holes describe are charac- 
terized by two energy scales, which we can choose to be yfp and T. The dimensionless ratio 
T/yfp parameterizes a family of symmetry-breaking black hole solutions in the $1-^0 plane. 
Let us consider p to be fixed. At some critical temperature, T = T c , the family of solutions 
intersects the flat line at tpQ = that corresponds to RNAdS. It is sometimes convenient 
to use T/T c instead of T/yfp to parameterize the broken solutions. In doing so, we should 
note that we hold p fixed: that is, what T/T c really means is (T / y/p) / (T c / y/p\) . This is 
distinct from what we would get by holding the chemical potential p fixed. Physically, we 
are studying the boundary theory at fixed charge density. 

Following [U E], we will also explore the conductivity for a selection of the solutions we 
generate. Because the condensate is purely s-wave, the conductivity is the same in any direc- 
tion: it is a complex scalar, a, whose real part measures the dissipative response and whose 
imaginary part measures the reactive response. One may calculate the conductivity using 
the formula <r{uS) = J(u)/E(u), where J is the current response to a spatially homogeneous 
electric field E. More precisely, if we introduce a perturbed gauge field, 

A^dx" = $dt + e~ iu;VJrot a x {r)dx , (17) 
and also perturb the metric by introducing an off-diagonal element 

gtx = e -^t e 2A(r )hUrh (lg) 

then the linearized equations of motion couple a x and ht x but do not require any additional 
fields to be perturbed. Allowing complex perturbations to the gauge field and the metric 
is a formal trick: it is understood that we take the real part at the end of the day. The 
linearized Einstein equations reduce to a single constraint between h tx and a x : 

h' tx + e- 2A &a x = 0. (19) 
This constraint can be used to eliminate h tx from the linearized Maxwell equations. The 
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resulting equation for A x is 



< + + - J < + - (^-J - 2 9 V --»)«. = <>. (20) 
Essentially this equation also appeared in [5]. The perturbation must obey standard 



horizon boundary conditions: it is infalling. With the requirements (12), this simply means 
that a x oc r~ lWv/ ^ for small r. At large r, a series solution in powers of e~ r suffices to show 
that 



a* r = a 



f+4V^ + .... (2i) 



The conductivity is 



^--^y, (22) 



simply because ai^e - ^^ 7 "* is the expectation value of the current J x in the dual field theory, 
while iua x °' } e~ luJyr ^° t is the electric field. 

Before delving into the results of numerically solving 11 ) , let us note an argument 



[2] in favor of the existence of symmetry-breaking black holes, for sufficiently large qL, 
that avoids numerics entirely and employs instead the Breitenlohner-Freedman (BF) bound 
[T2l IT3] j^] The starting point is to note that the T — > limit of the Reissner-Nordstrom black 
hole ^ has for its near- horizon geometry AdS-z x R 2 . Let's choose the radial variable r in 
(|2| so that the horizon is at r# = 0. Then at extremality one has 



(2tt) 3 L 3 r (2vr) 3 L 2 



P=7^7-^ ( 2 3) 



where e and p are defined as in (16). One easily sees that the radius of the AdS 2 geometry 
near r = is 

L*=±- (24) 

The transverse space is not significant for what follows (in the context of [2] it was S 2 not 
R 2 ). The point is to note that the BF bound points to an instability in AdS 2 when a scalar 
has a effective mass m 2 R satisfying 

m 2 m L 2 m <-\. (25) 



2 We note that a related line of argument has been refined and extended in [TU] to the case of neutral 
scalars with m 2 L 2 = —2. 
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Here m 2 R is the limit as r — > r# of the effective mass squared defined in [Tj: 

m 2 eS = m 2 + g tt q 2 $ 2 . (26) 



When the inequality (25 ) holds, according to the logic of [2], one should expect not just a few 
branches of symmetry-breaking solutions, but infinitely many, with arbitrarily many nodes 
in ip. If the inequality does not hold, it is not clear whether or not superconducting solutions 
exist: the question then depends on more detailed properties of the matter lagrangian. But 
at most one expects only finitely many branches of symmetry-breaking solutions. As usual, 
the branches with nodes can generally be expected to be unstable and/or thermodynamically 
disfavored. 

A short calculation based on (|2| and (26) shows that (25) is equivalent to 

m 2 L 2 -2q 2 L 2 <~. (27) 

For comparison, the BF bound in AdS^ is m 2 L 2 > —9/4. If m 2 L 2 < —9/4, then empty 
AdS± itself is unstable toward developing non-zero ip. If —9/4 < m 2 L 2 < —3/2, then 
AdSi is stable, but the argument explained in the previous paragraph indicates that there 
should be symmetry-breaking charged black hole solutions no matter what value q takes. If 
m 2 L 2 > —3/2, then only for sufficiently large q does this argument imply the existence of 
symmetry-breaking charged black hole solutions. Our numerics will focus on m 2 L 2 = —2 
and qL between 0.1 and 2^] 

3 Numerical results 



As explained following ( 13 ), the black hole solutions we are interested in can be parameterized 
by ipo (the value of the scalar field at the horizon) and $i (roughly, the electric field in the r 
direction at the horizon). The boundary condition on the scalar at the conformal boundary 
of anti-de Sitter space amounts to a single non-linear constraint between i/jq and $i. There is 
thus a one-parameter family of solutions. More precisely, there are several such families, or 
perhaps even infinitely many, with each family corresponding to black hole solutions where 



tp has zero, one, two, or more nodes. As explained following (14), we focus our attention on 
the family of solutions where ip is everywhere positive. To explore this family of solutions, 
we started by finding a few symmetry-breaking solutions with small These solutions 



3 We thank S. Hartnoll for correcting a numerical error in a previous version of this discussion. 
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sufficed to determine T c with good accuracy. (T c can also be determined by treating the 
scalar as a linearized perturbation of the Reissner-Nordstrom-anti-de Sitter solution, as in 
[I].) Then we proceeded to larger i/jq in small steps, ensuring at each step that the boundary 
condition on the scalar is satisfied. We iterated this process until we ran out of CPU time 
(the allocation was roughly 12 hours per branch) or until we encountered numerical errors 
suggesting that the solutions were no longer reliable. 

We further restrict our attention to the scalar potential (|3| with m 2 L 2 = —2. As ex- 
plained around ([6]), this corresponds to an operator in the dual field theory with dimen- 
sion 2, provided one makes the usual assignment of operator dimensions. There are three 
dimensionful parameters in the lagrangian: k, L, and q. But k enters only as an overall scale 
which doesn't affect the classical equations of motion or rescaled thermodynamic quantities. 
Only the dimensionless combination qL can enter into formulas describing the physics of the 
dual CFT, heuristically because conformal invariance prohibits any dimensionful scale in the 
theory. 

3.1 Thermodynamics 

Figure [T] shows a summary of thermodynamic quantities as a function of qL and the temper- 
ature. Of particular interest is Af, which is the rescaled free energy of the superconducting 
black hole, minus the rescaled free energy of the Reissner-Nordstrom black hole with the 
same temperature and charge density. The criterion for superconducting black holes to be 
preferred in the microcanonical ensemble is Af < 0. As is evident from figure [TJz>, Af is 
indeed negative for T < T c , no matter what qL is, over the range 0.1 < qL < 2. It appears 
that there are thermodynamically preferred superconducting black holes for all values of 
qL when the potential is ^ with m 2 L 2 = —2. This agrees with our BF bound argument 
from section [2} If we chose m 2 L 2 less negative, or positive, then we would expect there to 
be a minimum value of qL below which superconducting black holes either do not exist or 
are not thermodynamically favored. But, according to the arguments of [TJ, we expect that 
sufficiently large qL should produce superconducting black holes for any fixed value of m 2 L 2 . 

3.2 Index of refraction 

At least for qL ^ 0.7, a distinctive feature emerges at low temperatures: h has a "double- 
shelf" structure, visible in figures [2p and|2p. If the lower shelf extended infinitely far into the 
infrared, instead of terminating in a regular horizon, then we would wind up with a domain 
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Figure 1: (Color online) Thermodynamic quantities as functions of qL and temperature. In 
(A)-(D), the curves for each value of qL have different lengths because they terminate either 
where numerical errors began to cast doubt on the validity of the coldest solutions, or where 
we simply ran out of CPU time. All logarithms are natural logs. 
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Figure 2: (Color online) (A) The inverse index of refraction, l/rapR, as a function of qL. (B) 
The critical temperature T c and the inflection point temperature T* compared to the rescaled 
charge density p as functions of qL. (C) For qL = 0.7 and the coldest temperature we could 
reach, there is an inflection point in h(r) at the location of the black dot. In order to see a 
distinctive double shelf develop around this inflection point, we had to cool the black hole 
quite a bit more than we did for figure [l] (D) For qL = 1.4, there is a distinctive double-shelf 
structure in h(r) at a less extreme temperature than for qL = 0.7. The black dot is again at 
an inflection point of h(r). 

wall geometry similar to the one exhibited explicitly in |5iJ. In such a case, the infrared 
geometry would not be anti-de Sitter space, but rather some more general warped geometry 
with 50(2, 1) invariance. The overall geometry would then be a domain wall with no boost 
invariance separating two asymptotic geometries, both of which have SO(2, 1) invariance, 
but with different values of h, call them h\jy and /iir. The overall normalization of h can be 
changed by reparameterizing r, but the ratio 



r v 



h 



(28) 



IK 



is an invariant quantity. It can be termed an index of refraction, because it characterizes 
the maximum speed of transmission of signals in the infrared. Unlike in [5], we have been 
unable to explicitly construct zero-temperature solutions with an infinite shelf. However, 
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examination of a series of solutions with the same qL and decreasing values of T/T c shows 
that the lower shelf broadens out without changing its height appreciably. Because the 
existence of a shelf is a somewhat qualitative criterion, we studied instead the inflection 
points of h as a function of r. There are never more than two. If there are two, then the one 
at the smallest value of r is probably associated with double shelf behavior. We define rim 
for a finite-temperature solution by replacing hjR by h evaluated at its first inflection point. 
In figure |2]A, we show 1/niR for the coldest black hole we could construct as a function of 
qL. For qL ^ 0.7, we couldn't find solutions where there are two inflection points. For 
larger values, we observe that there is a temperature T* below which there are two inflection 
points, and above which there are not. We are not sure whether the double shelf appears at 
all values of qL or not. If it does, T* quickly becomes small, and tiir quickly becomes large, 
as qL decreases. Indeed, we found n IR « 4.8 x 10 5 at qL = 0.7. 



3.3 Conductivity 

For a selection of superconducting black holes, we computed the conductivity as a function of 



frequency by solving (20) with infalling boundary conditions at the horizon and using (22). 
Figure [3] shows the real and imaginary parts of the conductivity at various temperatures for 
qL = 0.7, 1.4, and 2. The results are roughly in line with the earlier results jl] in the probe 
approximation; also, they appear to approximately agree with the results of [10J, though the 
range of qL we surveyed is further from the probe approximation. 
Several qualitative points are worth noting: 

• The behavior Ima oc 1/u indicates the presence of a S(u) contribution to Re a. 

• When qL = 0.7, the conductivity curves appear to lie on top of each other. One way 
to understand this is that T c is very small compared to \fp for qL = 0.7, and making 
T even smaller doesn't significantly change the way the system responds. 

• The behavior of o at small uj/y/p is more "gap-like" for larger values of qL: that is, a 
is closer to for a wider interval. But we do not find evidence for a being strictly zero 
up to some finite u>. It could be that this happens for T exactly equal to 0; but we are 
inclined to think that there isn't a "true gap" in these systems. 

It is notable that a finite 5{uS) appears in Re ex even as T approaches T c — in fact, even above 
T c it persists. As noted in [TU], this is associated with the translation invariance of the 
system, which prevents a DC current from relaxing. Thus, infinite DC conductivity is not 
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Figure 3: (Color online) (A), (B) The real and imaginary parts of the conductivity for 
qL = 0.7 at various values of T/T c . (C), (D) The real and imaginary parts of the conductivity 
for qL = 1.4 at various values of T/T c . (E), (F) The real and imaginary parts of the 
conductivity for qL = 2 at various values of T/T c . 
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a feature that distinguishes the superconducting black holes from the normal state RNAdS 
black holes; rather, the breaking of the abelian gauge symmetry through the formation of 
the condensate is the distinguishing feature. 
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